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ABSTRACT 

We present the IR luminosity function derived from ultra-deep 70/im imaging of the GOODS-North 
field. The 70 /im observations are longward of the PAH and silicate features which complicate work in 
the MIR. We derive far- infrared luminosities for the 143 sources with Sjq > 2 mJy (S/N > Scr). The 
majority (81%) of the sources have spectroscopic redshifts, and photometric redshifts are calculated 
for the remainder. The IR luminosity function at four redshifts (z ~ 0.28, 0.48, 0.78, and 0.97) 
is derived and compared to the local one. There is considerable degeneracy between luminosity and 
density evolution. If the evolving luminosity function is described as p{L, z) ~ {l + zYp{L/ (1 -l-z)^, 0), 
we find q = — 2.19p-|-6.09. In the case of pure luminosity evolution, we find a best fit of p = 2.78j;Q32. 
This is consistent with the results from 24/im and 1.4GHz studies. Our results confirm the emerging 
picture of strong evolution in LIRGs and ULIRGs at 0.4 < z < 1.1, but we find no evidence of 
significant evolution in the sub-LIRG (L < IO^^Lq) population for z < 0.4. 
Subject headings: infrared: galaxies — galaxies: evolution 
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1. INTRODUCTION 

Deep mid-infrared surveys are revealing a population 
of mid and far-infrared luminous galaxies out to z ~ 3. 
These luminous (LIRGs, IO^^Lq < Ljr = is-ioooAim < 
IQI^Lq) and ultraluminous (ULIRGs, Lir > lO^^^^^) 
infrared galaxies are relatively rare in the local uni- 
verse, but become increasingly important at high red- 
shift, where dust enshrouded starb ursts dominate the to- 
tal cosmic star fo rmation rate (e.g. ICharv fc Elbad[2001l . 
iBlain et al.l[200l . 

The Infrared Space Observatory (ISO) showed that in- 
frared luminous starbursts were much more numerous at 
z ^ 1 than at the present time (jFranceschini et al]|2001l : 
lElbaz et al.l[200l . The ISO results were expanded upon 
by deep surveys at 24 fim with the Multiband Imag- 
ing P hotometer (MIPS) on the Spitzer Space T elescope 
(e.g. ICharv et al.l l2004l iPapovich et all l2004f ). Using 
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the excellent ancillary data in the Great Observatories 
Origins Deep Survey (GOODS) South and North fields, 
15 /im and total infrared luminosity fun ctions were de- 
rived from thousands of 24 /im so urces (|Le Floc'h et alj 
120051 iPerez-Gonzalez et al.l 120051 ). Strong evolution of 
the IR population was found and t he IR lumino sity func- 
tion evolves as (1 + ^)^ for z <l ()Le Floc'h e t al. 20051; 
iPerez-Gonzalez et al.|[2005l ). 

The 24 /im results are dependent on the set of SED 
templates used to extrapolate the 24 p,m flux densities 
to 15 fim and total infrared luminosities. Furthermore, 
significant variations in the bolometric correction are ex- 
pected as strong PAH and silicate emission and absorp- 
tion features are redshifted into the 24 /im band. Obser- 
vations with the 70 /im band of MIPS are closer to the 
peak in FIR emission and are not affected by PAH or 
silicate features for z < 5. They should therefore provide 
more robust estimates of the far-infrared (FIR) luminosi- 
ties. 

Studies by ISO in the FIR regime have been lim- 
ited in sensitivity (S'gopm ^ 100 mJy, Sirofim > 200 
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mJy) and redshift c o mpleteness (ISerieant et al.1 120041 : 
iTakeuchi et"all l2006f ). iFraver et all (|2006a[ ) derived a 
FIR luminosity function (LF) for the Extragalactic First 
Look Survey (xFLS) from Spitzer 70 /xm data, but 
this survey had incomplete redshift information at faint 
fluxes, and it was limited to z < 0.3 and bright (570^111 ^ 
50 mJy) sources. In this paper we present the infrared 
luminosity function up to redshift 1 from the ultra-deep 
70 fim survey of GOODS-N. 

We assume a Hubble constant of 71 kms^^Mpc^^, and 
a standard A-CDM cosmology with 17m = 0.27 and Ha = 
0.73 throughout this paper. We define the IR flux as the 
integrated fiux over the wavelength range 8 to 1000 /im. 

2. THE DATA 
2.1. Ultra- deep 70 fj,m Imaging 

The GOODS-N field is centered on the Hubble Deep 
Field North at 12h36m55s, -|-62°14ml5s. The MIPS 
70 ^m observations of GOODS-N we re carried out during 
Cycle 1 {Spitzer program ID 3325, (jFraver et al.ll2006bh 
and Cycle 3 (January 2006) for the Far Infrared Deep Ex- 
tragalactic Legacy project (FIDEL, Spitzer PID:30948, 
PI: Dickinson). Together these data map a region 10' x 
18' to a depth of 10.6 ksec. 

The raw data were processed off-line using the 
Germanium Reprocessing Tools (GeRT ) , follow ing the 
techniques described in iFraver et al.l ()2006bf) . We 
have cataloged 143 sources (over '^185 arcmin^) with 
S'7o>2.0mJy (S/N>3cr) in GOODS-N. The 70 ^im im- 
ages have a beam size of 18'.' 5 FWHM, and in the pres- 
ence of Gaussian noise the la positional error of sources 
is of the order °-^g™" , i.e. 3" for the faintest sources. 

2.2. Redshifts 

All 70 micron sources were matched to 24 micron and 
IRAC sources to obtain good positions. The best Spitzer 
position was then used to search for optical redshifts. 
About 7% of the 70 micron sources have more than one 
24 micron source within the 70 micron beam, and these 
were deblended individually (e.g. Huynh et al, ,200.7). 
Spectros copic redshifts are available for 116 o f the 143 
objects (|Cohen et al.l l2000t IWirth et all 120041 : Stern et 
al. in prep). 

Photometric redshifts were derived for 141 of the 
143 source s with th e exte nsive photometry available: 
ACS HST (iGiavalisco et aLll2004l) . U- (NOAO), BVRIz- 
(Subaru-SupremeCam) and JK- (NOAO/KittPeak- 
Flamingo) imaging. The photometric redshifts were cal- 
culated us ing the y 2 minimization technique as explained 
in Mobas her et al.l (2006). We have photometric red- 
shifts for 26/27 sources that don't have a spectroscopic 
redshift and we therefore have redshift information for 
142/143 sources. 

We quantified the reliability of the photometric red- 
shifts by examining the fractional error, A = (zphot — 
•Zspcc/(1 + -^spec)- For all 115 70 /im sources with both 
photometric and spectroscopic redshifts, we found the 
median fractional error. A, is 0.012 ± 0.20. Assuming 
the 6 cases where the fractional error is greater than 0.2 
are outliers, the success rate of the photometric redshift 
method is 95%. Removing the 6 outliers gives a median 
fractional error of 0.0014 ± 0.05. We therefore conclude 
that the photometric redshifts are statistically reliable. 



The 70 micron sources have a median redshift of 0.64 
(see Figure 1). The majority (79%) of sources lie at z < 
1, as expected for the survey sensitivity and steep k- 
correction that is present at 70 micron. 

3. INFRARED LUMINOSITIES 

Many authors argue that the MIR is a good indica- 
tor of the bolometric IR luminosity for normal and IR 
luminous galaxies (e.g. Chary and Elbaz 2001). Based 
on this, several authors have developed sets of galaxy 
templates that ca n be used to estimat e the total in- 
frare d luminosity (IChary fc Elbad 120011: iDale fc Heloul 
I200I iLagache et al.ll2003l i! 

We use the luminosity dependent SED templates based 
on local galaxies from Chary and Elbaz (2001) to deter- 
mine the IR luminosities of the 70 /im galaxies. However 
it is not clear whether local templates can accurately re- 
produce the MIR SED of distant galaxies because PAH 
and silicate absorption features are dependent on com- 
plex dust physics, including the intensity of the radiation 
field, the metallicity of the ISM, and the distribution of 
grain sizes. For this reason we determine the IR lumi- 
nosities of the 70 /im galaxies by fitting templates to the 
observed 70 /im flux density only, which is longward of 
the PAH and silicate features. 

The IR luminosities as a function of redshift are shown 
in Figure 1. Most of the sources below redshift z = 1 have 
LIRG-like luminosities. The higher redshift sources are 
luminous ULIRGs with possibly an embedded AGN. 

The estimated accuracy of the IR luminosity, from the 
70 /xm flux density calibration and PSF fltting errors 
alone, is 9%. However, the luminosities derived are de- 
pendent on the SEDs used. The adopted template SEDs 
do not reflect the full range of SEDs observed in galax- 
ies, and thus are the main source of systematic errors. 
For example, the total IR luminosity derived from the 
MIR regime can vary by a factor of 5 for local galaxies 
(|Dale et al.ll200l . We are working longward of the PAHs 
and silicate features which affected previous work based 
on the MIR, but, on the other hand, the restframe wave- 
lengths probed at 70 /im is affected by dust temperatures 
and emissivity. 

To test the consistency of our derived IR luminosi- 
ties and the application of the adopted SEDs, we use 
the well known FIR-radio correlation. The deep radio 
image of GOODS-N (5 /iJy rms at 1.4 GHz, Morrison 
et al. in preparation) detects 120/143 (84%) of the 70 
/im sources at 3ct or above. The FIR-radio correlation, 
q — log(FIR/S'i.4GHz), where 'FI R' here refers t o the 
flux between 40 and 120 /im (e.g.lYun et al."2001'), has 
an observed local value oiq = 2.34±0.3 (Yun et al. 2001). 
Adop ting an average fact or of 2.0 between IR and FIR 
(e.g. iDale fc Heloul l200^ and a radio spectral index of 
a = —0.8^, we find g = 2.2 ± 0.2 for the radio detected 
sources. Including the 24 /im data in the fits to the SEDs 
gave a slightly larger dispersion in q. This suggests that 
the IR luminosities as estimated from the 70 /im data 
alone are reasonable. 

4. INFRARED LUMINOSITY FUNCTIONS 

In this Section we explore the evolution of the IR lu- 
minosity function between redshifts and 1. 

2 oc y-^. 
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4.1. Methodology 

The luminosity functions were derived for 4 redshift 
bins, 0.2 < z < 0.4, 0.4 < z < 0.6, 0.6 < z < 0.9, and 
0.9 < z < 1.1 using the usual 1/Knax method (Schmidt 
1968). These redshift bins were made wide enough so 
that there is a reasonable number of sources for calculat- 
ing the luminosity function. The bins have median red- 
shifts of 0.28, 0.48, 0.78, and 0.97, so a moderate range 
in redshift is explored. The comoving volume for each 
source is Vmax — yzrc^^y, — ^z„i„, where Zmin is the lower 
limit of the redshift bin, and Zmax is the maximum red- 
shift at which the source would be included in the cata- 
log, given the limiting 3ct limit, or the maximum redshift 
of the bin. 

As mentioned in Section 2.2, we have almost complete 
redshift information on the 70 /im sample. A correction 
factor for each individual source was computed to correct 
for source detectability across the full image and flux 
boosting (i.e. the over- estimation of the flux densities of 
low SN sources) . This correction was calculated us i ng the 
Monte Carlo approach described by iCharv et al.l ()2004f ) 
and it is the same co rrection applied to the source counts 
(|Fraver et al.ll2006H) . 

4.2. Results and Discussion 

The luminosity functions were derived from the rest- 
frame IR /im luminosities. In Figure 2 we plot the lumi- 
nosity functions for the redshift bins explored, and the 
data is summarized in Table 1. The local IR lumin os- 
ity function from IRAS sources () Sanders et al. I [20031 ) is 
plotted for comparison. 

For each luminosity bin the uncertainties, Up, were es- 
timated using the Poisson statistics on the number of 

sources, so dp = ij^ ) ■ Monte Carlo simulations 

were also performed to disentangle the uncertainties in 
the derivation of the luminosity function due to photo- 
metric errors. Each source was randomly given an IR 
luminosity within the uncertainty estimates and the lu- 
minosity function was re- calculated. We find this adds 
between 0.03 to 0.09 dex to the luminosity function un- 
certainty, depending on the bin, but the Poisson statistics 
dominate the uncertainties. 

We do not find any significant evolution in the sub- 
LIRG population {L < IO^^Lq) for the lowest redshift 
bin {z < 0.4) (Figure 2). The high redshift LFs show 
evidence for strong evolution of LIRGs and ULIRGs at 
z > 0.4. 

The IR LFs derived here are co nsistent with that 
derived from 24 micron ()Le Floc'h e t al. 2005) for the 
overlapping luminosity bins at z < 1. This implies 
that, on average, similar bolometric luminosities are de- 
rived from 24 and 70 micron for moderate luminosity 
{L < and moderate redshift sources {z < 1) 

sources. We can not say if this i s the case for high lu- 
minosities and h igh redshifts (e.g. IChapman et al.lT2005t 
iPope et al]l2006) . as those sources are rare in the 70 mi- 
cron data. Recent 70 micron stacking analysis of galaxies 
at z ~ 2 show that 24 micron observations at high red- 
shift over-estimate LI R in comparison to 70 micron and 
other LIR indicators (jPaddi et al.ll2007t IPapovich et al.l 
l2007l) . 

To explore the evolution of IR sources we use the an- 



alytical form of the local luminosity function (LF) from 
ISanders et al.l (|2003( ) that comprises of a double power 
law: p cx for log(L/L0) < 10.5, p (x L'^ ^ 

for log(L/_L0) > 10.5. We assume that the evolv- 
ing luminosity function can be described by p{L, z) = 
g{z)p[L/ f{z),Q\. In this sense g{z) and f[z) describe 
the density and luminosity evolution of the LF, respec- 
tively. The commonly used form of evoluti on is to assume 
f(z) = (l + z)P and q(z) = (1 + z)« (e.g. lCondonl[l98l 
iHaarsma eraLll2000[ ). 

Using minimization, we examine the best fit to the 
evolution of the IR LF. There is a well known degeneracy 
between density and luminosity evolution. We find that 
the best fit evolution parameters follow the relation q — 
— 2.19p+ 6.09. In the case of pure luminosity evolution 
{q = 0), we find p = 2.78toit- 

These evolution constraints are broadly consistent with 
24 iim studies which found p = 3.2^0^2 and q = 0.7'!i'^'^ 
for t he infrared luminosity function (|Le Floc'h et al.l 
l2005f ). Our results are also in good agree ment with previ- 
o us studies of IR sources (e.g. iFranceschini et al.|[2b01i ). 

iHopkinsI ()2004f ) combined star formation rate data with 
faint radio source counts to find p = 2.7 ± 0.6 and q = 
0.15 ± 0.60. If only pure luminosity ev olution of radio 
sourc es is conside red then p = 2.5 ib 0.5 (jSevmour et al.l 
120041 ) or p = 2.7 (jHuvnh et al.ll2005f ). So our results are 
consistent with constraints on the evolution of the star 
forming population from deep radio surveys, indicating 
that the radio sources overlap with the ultra-deep 70 fim 
population, as expected. 

The constraints on the evolution of the IR LF can be 
used to determine the cosmic star formation r ate (S FR) 
density. Using the calibration from lKennicuttI (jl998l ) and 
integrating over galaxies with 8.5 < log(L/L0) < 12.5, 
we find the SFR density at z = 1 is 0.15t5J;J^3 Mq yr"! 
Mpc~'^ for the best fit pure luminosity evolution case. 
Here the uncertainties in SFR density do not include 
the systematics in the FIR/SFR calibration, which add 
about 0.3 dex to the absolute uncertainty. The SFR 
density derived here is lower than that e stima ted by the 
evolutionary models of ICharv fc Elba3 (j2001h by about 
a factor of 1.7, but it is con sistent with extinct ion cor- 
rected optic al measures (e.g iKewlev et al.ll2004[) and 24 
pm resuhs (Le Flo c'h et al.ll2005f ). 

The AGN in our sample can be identified us- 
ing the__deeg__2_Ms JC^ay observations of GOODS- 
N (j Alexander et al.l l2003f l. Sources are classed as X- 
Ray AGN from X-ray band ratios . X-ray luminosity , 
and X-ray-to-opt ical ilux ratios ( Alexander et al.l l2003l : 
iBauer et al1l2004f ). At redshifts z < 0.6 we find only 7% 
of the 70 ^m sources are X-Ray AGN, but this fraction 
increases to 27% for the 0.9 < z < 1.1 redshift bin. The 
highest redshift LF in Figure 2 is contaminated by X-Ray 
AGN but this does not significantly affect the evolution 
derived in this work. 

5. CONCLUDING REMARKS 

Based on ultra-deep 70 /xm observations of GOODS- 
N, and the spectroscopic and photometric redshifts avail- 
able of galaxies in this well-studied field, we have derived 
luminosity functions for z = 0.3 to z = 1.1. We find 
strong evolution in galaxies with Lm > 10^^ Lq at red- 
shifts z > 0.4. Assuming pure luminosity evolution of 
the form (1 + z)p, we find p = 2.78toit This confirms 
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TABLE 1 

The IR luminosity functions derived from the 1/Vm 



ANALYSIS. 
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9 1.24 ± 0.41 X 10-3 
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0.6 < 2 < 0.9 




11.10 


11.40 


11.25 


9 1.17 ± 0.47 X 10-3 


0.05 
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11.47 


13 7.90 ± 2.31 X 10-4 
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12.00 
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7 2.48 ± 0.94 X lO-'' 
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8 1.11 ± 0.55 X 10-^ 


0.06 


11.70 
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"The errors quoted arc Poissoii uncertainties. 

''This is the additional uncertainty to be added to the LF, calculated from Monte Carlo simulations of the uncertainty in the IR luminosity 
of each source. See Section 4.2. 



the strong evolution in LIRGs and ULIRGs between red- 
shift and 1 that has been seen in previous work. The 
depth of the 70 /ini data allows us to probe sub-LIRG lu- 
minosities, and we find little evolution in this population 
for z < 0.4. In the case of pure luminosity evolution, we 
find the star formation rate density at = 1 is 0.151q q3 
Mq yr-i Mpc"^. 

This is the first result from an ultra-deep FIR survey 
that reaches ^; ^ 1. However, we are limited by poor 
statistics - the number of bins available for the LF at 
each redshift is limited by the small number of cataloged 
sources. The area covered is only 10' x 18' so these results 
are also affected by cosmic variance. 

The Far Infrared Deep Extragalactic Legacy (FIDEL) 
Spitzer legacy project, currently underway, will cover the 
extended Chandra Deep Field South and the extended 



Groth Strip at 70 /im with similar depths to the GOODS- 
N. The total area covered will be about 9 times that used 
in this work, and the FIDEL project will detect over 1000 
LIRGs at moderate redshift. So in the near future, large 
ultra-deep FIR surveys such as FIDEL will enable even 
more detailed studies of the FIR luminosity function and 
the evolution of infrared galaxies. 
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was provided by NASA through an award issued by 
JPL/Caltech. 



REFERENCES 



Alexander, D. M., et al. 2003, AJ, 126, 539 

Bauer, F. E., Alexander, D. M., Brandt, W. N., Schneider, D. P., 
Treister, E., Hornschemeier, A. E., & Garmire, G. R 2004, AJ, 
128, 2048 

Blain, A. W., Smail, I., Ivison, R. J., Kneib, J.-P., & Frayer, D. T. 

2002, Phys. Rep., 369, 111 
Chapman, S. C, Blain, A. W., Smail, I., & Ivison, R. J. 2005, ApJ, 

622, 772 

Chary, R., et al. 2004, ApJS, 154, 80 
Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 

Cohen, J. G., Hogg, D. W., Blandford, R., Cowie, L. L., Hu, E., 

Songaila, A., Shopbell, P., & Richberg, K. 2000, ApJ, 538, 29 
Condon, J. J. 1984, ApJ, 284, 44 
Daddi, E., et al. 2007, ArXiv c-prints, 705 
Dale, D. A., et al. 2005, ApJ, 633, 857 
Dale, D. A., & Helou, G. 2002, ApJ, 576, 159 

Elbaz, D., Ccsarsky, C. J., Chanial, P., Aussel, H., Franceschini, 

A., Fadda, D., & Chary, R. R. 2002, A&A, 384, 848 
Franceschini, A., Aussel, H., Ccsarsky, C. J., Elbaz, D., &; Fadda, 

D. 2001, A&A, 378, 1 

Frayer, D. T., et al. 2006a, AJ, 131, 250 
Frayer, D. T., et al. 2006b, ApJ, 647, L9 
Giavalisco, M., et al. 2004, ApJ, 600, L93 

Haarsma, D. B., Partridge, R. B., Windhorst, R. A., & Richards, 

E. A. 2000, ApJ, 544, 641 



Hopkins, A. M. 2004, ApJ, 615, 209 

Huynh, M. T., Jackson, G. A., Norris, R. P., & Prandoni, I. 2005, 
AJ, 130, 1373 

Huynh, M. T., Pope, A., Prayer, D. T., & Scott, D. 2007, ApJ, 659, 
305 

Kennicutt, R. C, Jr. 1998, ARA&A, 36, 189 

Kcwlcy, L. J., GcUer, M. J., & Jansen, R. A. 2004, AJ, 127, 2002 

Lagachc, G., Dole, H., & Puget, J.-L. 2003, MNRAS, 338, 555 

Lc Floc'h, E., et al. 2005, ApJ, 632, 169 

Mobasher, B., et al. 2006, ArXiv Astrophysics e-prints 

Papovich, C., et al. 2004, ApJS, 154, 70 

Papovich, C, et al. 2007, ArXiv e-prints, 706 

Perez-Gonzalez, P. G., et al. 2005, ApJ, 630, 82 

Pope, A., et al. 2006, MNRAS, 370, 1185 

Sanders, D. B., MazzarcUa, J. M., Kim, D.-C., Surace, J. A., & 

Soifer, B. T. 2003, AJ, 126, 1607 
Serjeant, S., et al. 2004, MNRAS, 355, 813 

Seymour, N., McHardy, I. M., & Gunn, K. F. 2004, MNRAS, 352, 
131 

Takcuchi, T. T., Ishii, T. T., Dole, H., Dennefeld, M., Lagache, G., 

& Puget, J.-L. 2006, A&A, 448, 525 
Wirth, G. D., et al. 2004, AJ, 127, 3121 

Yun, M. S., Reddy, N. A., & Condon, J. J. 2001, ApJ, 554, 803 



5 




Log (L„) [Lg] 

Fig. 2. — The IR luminosity function (LF). Crosses mark the local LF from Sanders et al. (2003) and the corresponding solid line is the 
double power fit to the local data. The symbols mark the LF calculated in this work at redshift 0.3 (upside down triangles), 0.5 (squares), 
0.8 (circles) and 1.0 (triangles). The lines are the local LF evolved to the corresponding redshift with the best fit pure evolution parameters. 
The arrows indicate bins which are incomplete because of the survey sensitivity. The horizontal error bars indicate the binsizes. 



